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Abstract Large-scale bulk motions and hydrodynamic turbulence in the intergalactic 
gas inside clusters of galaxies significantly broaden X-ray emission lines. For lines of heavy 
ions (primarily helium-like and hydrogen-like iron ions), the hydrodynamic broadening is 
significantly larger than the thermal broadening. Since cluster of galaxies have a negligible 
optical depth for resonant scattering in forbidden and intercombination lines of these 
■ ions, these lines are not additionally broadened. At the same time, they are very intense, 

which allows deviations of the spectrum from the Gaussian spectrum in the line wings to 
be investigated. The line shape becomes an important indicator of bulk hydrodynamic 
processes because the cryogenic detectors of new generation of X-ray observatories will 
^ ■ have a high energy resolution (from 5 eV for ASTRO-E2 to 1-2 eV for Constellation-X 

and XEUS). We use the spectral representation of a Kolmogorov cascade in the inertial 
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range to calculate the characteristic shapes of X-ray lines. Significant deviations in the 



line profiles from the Gaussian profile (shape asymmetry, additional peaks, sharp breaks 



in the exponential tails) are expected for large-scale turbulence. The kinematic SZ effect 



and the X-ray line profile carry different information about the hydrodynamic velocity 
distribution in clusters of galaxies and complement each other, allowing the redshift, the 
peculiar velocity of the cluster, and the bulk velocity dispersion to be measured and 
separated. 
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1 The General Picture 

The hot intergalactic gas in clusters of galaxies forms an extended atmosphere in the 
gravitational potential well produced mainly by the weakly interacting dark matter, whose 
mass exceeds the intergalactic gas mass by a factor of approximately six. The cluster gas 
temperature reaches 2-10 KeV. The speed of sound in this gas is 1000-1500 km s _1 . 
Galaxies move through the cluster gas with subsonic, sonic and supersonic velocities. 
Energetically, the presence of galaxies is of little importance, because the total mass 
of the galaxies is appreciably smaller than the total cluster mass. The intergalactic gas 
consists of completely ionized hydrogen together with ~ 25 % (by mass) of helium. Heavy 
elements, including iron, are represented at a 30 - 40% level of their solar abundance. 

Due to the growth of large-scale cosmological density perturbations, clusters of galax- 
ies sometimes merge together and capture surrounding galaxies and groups of galaxies. 
During such merger, the colliding components move with supersonic or transonic veloci- 
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ties. Each merger is accompanied by shock wave formation and turbulence generated on 
smaller scales. The characteristic time between mergers is long. The effects of no more 
than one to three large mergers appear to be observed simultaneously; a large merger is 
a merging of clusters with comparable masses. The capture of low-mass companions or 
individual galaxies gives a smaller contribution to the generation of turbulence. 

It follows from iV-body simulations with cold dark matter and gas that the turbulent 
pressure can account for up to 15% of the thermal pressure of the intergalactic gas in 
relaxed clusters of galaxies. This value implies very high turbulent pulsation velocities 
(up to 40% of the speed of sound). Remarkably, the calculations by different groups yield 
similar results (the methods, codes, and results were carefully compared by Frenk et al. 
1999; see also Norman and Bryan 1999). 

Activity of accreting supermassive black holes in the dominant galaxy inside the cluster 
can trigger additional small-scale turbulence due to existence of collimated high-velocity 
mass outflows and interaction of intergalactic gas with relativistic jet (Churazov et al. 
2002a, 2002b). 



1.1 X-ray Line Profiles 

In 2005, IS AS and NASA are planning to launch the ASTRO-E2 satellite 1 . This satellite 
will be equipped with cryogenically cooled X-ray bolometers. They will be placed in the 
focal planes of grazing-incidence telescopes. The expected energy resolution will be ~ 
5 eV (FWHM) in the photon energy range from 0.5 to 10 keV, which includes the line 
of the helium-like iron ion (the iron ion with two electrons) at energy 6.7 keV. Since 
below we compare the line profiles with Gaussian fits, it is important to note (Porter and 
Mitsuda 2003) that this resolution corresponds to a dispersion of 2.7 eV. The prospective 
Constellation- X 2 and XEUS 3 missions will have X-ray detectors with an energy resolution 
as high as 1-2 eV. Therefore, it will become possible to detect line shifts and broadenings 
in rich cluster of galaxies that correspond to turbulent velocities ~ 3 x 10~ 4 of the speed 
of light, i.e. up to 100 km s -1 or less than 10% of the speed of sound. 

During turbulent pulsations, iron ions move in a mixture with hydrogen and helium 
nuclei. Accordingly, all ions have the same hydrodynamic velocities. The thermal veloci- 
ties of ions of different types with the same temperature greatly differ because of the large 
difference between the nuclear masses. Therefore, the Doppler turbulent broadening of 
iron lines AE/E ~ v tUT h/c can significantly exceed their thermal broadening. The thermal 
velocity of iron ions is ~ ^Jm v jmF e = 1/ V^Q = 13% of the speed of sound, where m p and 
mpe are the masses of the proton and the iron atomic nuclei, respectively. The amplitude 
of the turbulent velocity pulsations produced by a cluster merger can exceed the thermal 
velocities of iron ions several fold. Thus, iron ions become effective tracers of turbulent 
velocity fields in the cluster interiors. 

Let us introduce the parameter 



Ma, = = V2p^ = A.2 rr~~~~T~\ \ ^ (^) 



v /fcTA^ Az/ dt • 300 [km s- 1 ] V kT \ 56m p 



1 http : / / www . isas . ac . j p/e /enterp / missions / ast ro-e2 / 
5 http://constellation.gsfc.nasa.gov/docs/main.html 
"http://astro.estec.esa.nl/SA-general/Projects/XEUS/ 
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that characterizes the ratio of the turbulent (iWb ~ u rms ) and thermal velocity scales for 
an ion of mass m^. It resembles the Mach number: Ma = u rms / c s (the ionic Mach number), 
where c s is the speed of sound. Parameter (1.1) defines the ratio of the hydrodynamic 
(Az/ DH ) and thermal (Az/ DT ) Doppler broadenings: 

A «rm S a ^2kT/rrn 

Az/ DH = , Az/ DT = vo • (1.2) 

c c 

In (1.1) and (1.2) u rms is the line-of sight or radial hydrodynamic velocity dispersion and 
u Q is the frequency at the center of the line profile. 

The most intense lines observed in the spectral range between 2 and 10 keV are those 
of helium and hydrogen-like iron ions. The equivalent width of the most intense iron lines 
reaches 100-500 eV. At the thermal Doppler width of ~ 3 eV, they must rise by tens and 
hundreds of times above the smooth continuum associated with hydrogen-helium plasma 
bremsstrahlung. XMM observations of the central part of the Perseus cluster of galaxies 
gave 10 5 photons in 50 ks in the complex of iron K-lines near 6.7 keV (Churazov et al. 
2003). The Constellation-X and XEUS satellites will have an effective area that is tens of 
times larger than that of XMM, which will make it possible to study in detail the weak 
wings of X-ray lines. 

The spectral surface brightness of a cluster in a line is given by the formula 



I(AE,y,z) = J dxj{E ,T e (x,y,z)}n e (x,y,z) 



drii(x,y,z, AE) 
d(AE) 



where / is measured in photons cm -2 s _1 eV -1 (square angular minute) -1 , the emissivity 
coefficient j characterizes the rate of production of excited ions and emission of photons 
with energy E by these ions, n e is the electron density, and drtij 'd(AE) is the number 
of ions of a given type in the range of velocities that shift the line photons to the energy 
range from AEtoAE + d(AE) (in the frame of reference associated with the ion, the pho- 
ton energy is Eq, without Doppler and intrinsic broadenings). The energy AE = E — Eq 
is measured from the line center and is proportional (see below) to the line-of-sight ion 
velocity. The line-of-sight velocity of ions can be determined by taking into account their 
thermal and hydrodynamic velocities. We assume that the thermal velocity distribution 
of ions is Maxwellian with an ion temperature equal to the electron temperature. The 
hydrodynamic velocity distribution can be calculated in terms of the turbulence model 
presented below. In the above formula, the natural line width is disregarded (see for- 
mulas (12.2) - (12.5) in the full paper). Emissivity coefficient j depends on the electron 
temperature T e . In an isothermal plasma at a constant (in volume) elemental abundance, 



n e rii oc n 2 e . 



In this paper, we consider qualitative effects, primarily the line broadening, its spec- 
tral shift due to bulk and turbulent motions and the degree of deviation of the emission 
line profile from the Gaussian profile. Therefore, for illustration, we investigate the sim- 
plest case of an isothermal cloud with constant density and elemental abundance over its 
volume. In this case, the above formula for the brightness / takes the form 

I(AE,y,z) = j{E ,T e )n e J dx ~^^y = a 3 n t J dxip(x,y,z, AE), 
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where a is the elemental abundance, and the function ip specifies the local line-of-sight 
velocity distribution. We ignore the nonlinear density variations due to turbulent pulsa- 
tions. 

An efficient method for solving the problem of the iron line profile involves direct 
numerical simulation. The first attempt of this kind was made by Sunyaev, Norman, and, 
Bryan (2003). Below, we attempt to construct a relatively simple turbulence model (see 
formula (4.1) below) to understand in which cases the shape of a turbulently broadened 
line can differ from the Gaussian shape expected from the central limit theorem of the 
probability theory. The model is based on the assumption about a Kolmogorov cascade 
in the inertial range limited by the mixer size L m j x and the viscous scale L v . The number 
of mixers in the large-scale case (-Z^mix ~ L) is limited. As a result, the line spectrum 
significantly deviates from the Gaussian spectrum 4 . 

Fig. shows the characteristic profile of the helium-like iron (Fe XXV) permitted 
w— line (curve 4) obtained in the adopted model. As we see (cf. curves 2 and 4), there 
are qualitative deviations from the commonly assumed Gaussian line profile 

Mrms exp [-(1/2) u 2 /(u 2 rms + fcT/mQ] _ (2n)-^ f (1/2) u 2 \ 

with a dispersion the square of which is equal to the sum of the squares of the thermal 
and hydrodynamic dispersions (curve 2). In what follows 

u = u/u rms , AE[eV] = -6.7 [eV] (u^/SOOfkms -1 ]) u, (1.3) 

u is a dimensionless velocity (in units of dispersion), and AE is the deviation from the 
profile center. Formula (1.3) is used to recalculate u to AE. Parabola 1 in Fig cor- 
responds to the Maxwellian distribution of ions f T {u) = (Ma,/V^) exp(- Ma? u 2 /2) 
(purely thermal broadening). All distributions where normalized to unity j_ f(u)du = 
1. Curves 3, /ht(w), and 4, /htl(w), represent the turbulent line profiles. The distribution 
/ht(w) includes hydrodynamic (H) and thermal (T) broadenings, while the distribution 
/htl(w) includes hydrodynamic, thermal, and Lorentz (L) broadenings. The intrinsic or 
Lorentz broadening was calculated for the Fe XXV w— line 5 . 

Fig. |2]shows, first, the computed spectra (fluctuating solid curves), second, the thermal 
Gaussians of the helium-like iron w-line (dashed curves, narrow parabolas), and, third, 
the Gaussian fits to the computed spectra (longOdashed wide parabolas). These spectra 
were constructed for nine points of the cluster projection onto the plane of the sky. The 
Gaussian fits are specified by two parameters (shift U\ and broadening mht) that can be 
determined from the computed spectrum (see Table 1). These nine columns correspond 
to the center, vertices, and middles of the sides of the square in the plane of the sky 
described in Sections 1.2 and 1.3, see also Section 11 of the full paper. Note that curve 
3 in Fig. ^corresponds to the line profile shown in the middle panel in the upper row of 

Fig. m 

4 The influence of the number of mixers on the deviation from Gaussian distribution is discussed in 
Section 7, see also the full paper. Large-scale motion shifts the line center, Fig. 2 and Table 1 below. 
The presence of a Kolmogorov tail in small scales does not imply that the profile will be Gaussian (see 
Table 2). 

5 See calculations in Sections 12 and 13 in the full paper. 
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Figure 1: Fe XXV w-line spectrum (4) at kT = 3 keV; 1 - the thermally broadened 
line profile (without hydrodynamic broadening), 2 - the Gaussian profile with a total 
(hydrodynamic plus thermal) dispersion mht = 190 km s -1 , and 3 - the line profile 
without Lorentz broadening. 



Fig. El (a) shows the radiation spectrum for a hot plasma with a temperature of 3 keV 
and a normal chemical abundance 6 in the energy range 7 near the permitted helium-like 
iron line at w 6.7 keV. The complex of iron lines near the hydrogen-like iron line (E w 
6.9 keV) is shown in Fig. El (b) at kT = 8 keV. These plots show the individual lines 8 in 
the form of narrow peaks with a width close to 10 -2 eV. The plots were constructed by 
using the APEC code (Astrophysical Plasma Emission Code; Smith et al. 2001), which is 
part of the XSPEC VI 1.2 code. All lines were normalized to the most intense permitted 
w-line (Eq = 6.7005 keV). As an illustration, the dashed curves in Fig. 01 (a) and (b) 
indicate the same spectra broadened by thermal motions of the ions. 

Fig. El (a)-(c) show the complex of lines near 6.7 keV taking into account both ther- 
mal and hydrodynamic motions. These spectra were computed by folding the profiles 
presented in the upper middle, central, and lower left panels of Fig. [2] with the complex 
of individual lines and satellites shown in Fig. El (a). The dashed curves again repre- 
sent the broadening of the individual lines due to thermal ion motions in the absence of 
hydrodynamic motions. 

The solid and dashed curves in Figs. 0] (a)-(c) were normalized to the same number 
of photons in the whole complex of lines. 

The plasma emission is a function of the temperature. An increase in temperature 
affects both the relative line intensity and the thermal line width. At high temperatures 
(e.g., kT = 8 keV), the complex of hydrogen-like iron lines near energy 6.9 keV becomes 
intense, Fig. El (b). 

Using Figs. El (a) and (b), we wish to emphasize the following. The most favorable 
energy range for observing hydrodynamic iron line broadening effects is to the right of the 
w-line center, because there are virtually no other intense lines at kT = 3-10 keV within 
several tens of eV of the center of this line. Thus, the far wings of the broadened line can 
be studied in detail. 

Of considerable interest is also the iron z-line at photon energy Eq = 6.6366 keV that 
corresponds to the transition from the ls2s triplet to the ground state: 3 S — s- 1 S. In 
this case, the left line wing is favorable for observations. It is possible to simulate the 
spectrum of the entire set of lines in both complexes and to compare it with the observed 
spectrum. 

Fig. 0] (c) shows a very large shift of the line observed from the column in the lower left 
panel in Fig. El All of this column flies toward us. The velocity dispersion in it slightly 
exceeds the thermal velocity dispersion. In contrast, the upper right panel of Fig. E] shows 
that the corresponding column flies away from us with a slightly lower radial velocity and 
larger velocity dispersion. Collectively, these two panels suggest the presence of a large- 
scale motion similar to the rotation of the entire cluster gas - this is a manifestation of a 
large-scale mixer, i.e., the last large merger. 

In the central part of the cluster (Figs. El and 0] (b)), the broadening and the shift 
are so large that the individual lines and satellites in the complex of lines are barely 
discernible. 



6 The abundance influences only the relative continuum intensity. 

7 The helium-like iron lines and their satellites are intense in this energy range. 

8 Without including their radiation widths. 
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Figure 2: Helium-like iron line profiles at nine points in the image of a cluster with 
kT = 3 KeV. The solid curves represent line profiles shifted and broadened due to hy- 
drodynamic and thermal motions and the dashed curves (narrow parabolas) specify the 
thermal line broadening. The hydrodynamic and thermal velocity dispersions and the 
velocity corresponding to the shift of the line profile centroid are given in Table 1 for each 
of the nine lines. The long-dashed curves (wide parabolas) represent the Gaussian fits to 
the computed spectra. 



7 



T=3 keV 



T=8 keV 




6.61 6.62 6.63 6.64 6.65 6.66 6.67 6.68 6.69 6.7 6.71 6.91 6.92 6.93 6.94 6.95 6.96 6.97 6.98 6.99 

Energy (keV) Energy (keV) 



Figure 3: Spectra for the complexes of iron lines near 6.7 keV at kT e = 3 keV (a) and 
near 6.9 keV at kT e = 8keV (b). The dashed curves specify the thermal broadening of the 
lines. The lines were normalized to the most intense permitted helium-like iron w-line. 

Table 1. Shifts and broadening of the X-ray line profile at Eq = 6.7 keV. The positions of the 
cells in the table and the data in them correspond to the spectra of the nine line-of-sight columns 
presented in Fig. El The fourth row of the table gives the velocity dispersion for the cluster as 
a whole. This is dispersion of the integrated spectrum in Fig. \5\ In the Table, u\ is the shift of 
the line profile center in [km s _1 ], and uht is the dispersion of the computed distributions from 
Fig. El 
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Since the ASTRO-E2 spectrometers have a limited angular resolution, they will yield 
detailed images only for the nearest rich clusters of galaxies. For distant clusters, the line 
profiles from the entire cluster will be investigated. Figs. 0] (d) and 03 give the first idea 
of this profile. A comparison of Fig. 0] (d) with Figs. 0] (a)-(c) shows how informative 
the set of high resolution spectroscopic data is even with a limited angular resolution. 
The profile in Fig. El was obtained by averaging the nine profiles presented in Fig. El with 
equal statistical weights. Recall that each of these nine profiles was obtained for a very 
narrow unit column in a cube with a volume of 160 3 mesh points. Clearly, the line profile 
from the entire cube must be even smoother. On the other hand, the shifts of the profiles 
shown in Fig. Eland Table 1 suggest the presence of intense large-scale motions inside the 
cluster. In observations of the cluster as a whole, they give a large contribution to the 
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Figure 4: (a,b,c) Spectra for the complex of iron lines near 6.7 keV for kT e = 3 keV are 
obtained by folding with the profiles shown in the upper middle, central, and lower left 
panels of Fig. |2l respectively. The dashed curves specify the spectrum for the thermally 
broadened lines of the complex. All curves were normalized to the same total number of 
photons in the complex of lines, (d) An approximate spectrum for the complex of iron 
lines near 6.7 keV from the cluster as a whole obtained by folding the profile in Fig. |3] 
with the spectrum in Fig. |3] (a). 
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dispersion but must leave traces in high-energy-resolution spectra. The profile in Fig. El 
is much broader than any of the nine profiles in Fig. |21 (see Table 1); in particular, it is 
broader than the profile in Fig. ^ This is the result of large-scale motions. 



r Photons/(cm 2 s eV) 




AE, eV 

Figure 5: The iron line profile (curve 3) at Eq = 6.7 keV obtained by summing the 
nine profiles shown in Fig. 121 This profile was also normalized to the total number of 
photons. To a first approximation, this profile corresponds to the line profile observed 
from the entire cluster. The thermal broadening (curve 1) and the two Gaussian curves 
for the total broadening in the entire cluster (total dispersion 460 km s _1 , curve 4) and 
for dispersion w rms = 270 km s _1 (curve 2) corresponding to the central panel in Fig. |2] 
are shown for comparison. 

The spectrum in Fig. 0] (d) was obtained by folding the profile in Fig. El with the 
complex of iron lines near 6.7 keV shown in Fig. El (a). We see individual spectral features, 
but their amplitude is appreciably smaller than that predicted in the model with thermal 
broadening and a smaller hydrodynamic velocity dispersion. The high-energy wing of the 
permitted helium-like iron w-line makes it possible to judge the total velocity dispersion 
in the cluster. 

The turbulent broadening determined for all the nine profiles shown in Fig. |21 cor- 
responds to a large (at least severalfold) decrease of the optical depth in the resonance 
w-line. That is why we disregard the broadening of this line due to resonant scattering 
(Gilfanov et al. 1987) in the plots of Fig. 2-5. 

To elucidate the turbulence model, Fig El shows a Kolmogorov spectrum of the Fourier 
velocity amplitudes u n with a power-law scaling in the wavelength range L v < A < L m i x < 
L and harmonic numbers 1 < iV m i x < n < N v , whereL is the cluster scale, X = X n = L/n, 
k = k n = 2n/\ n = 2Trn/L, N mix = L/L mix , N v = L/L u , and L u is the viscous damping 
scale. Note that the total number of energy-containing eddies in volume L 3 is iV^ ix . 
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ln(X/L) = ln(l/«) 




Figure 6: The A and k fluctuation spectrum in a Kolmogorov model. The scales {L,l), 
(L mix , iV mix ), and (L U ,N U ) characterize the cluster, the turbulence-maintaining mixers, 
and the viscous dissipation region, respectively. 

Below, we show that the velocity dispersion and the hydrodynamic Doppler broadening 
(1.2) are determined by large-scale fluctuations with A ~ L m i x . The fractal line profile 
fluctuations, which are smoothed out by thermal broadening, are associated with the 
Kolmogorov small-scale tail. 

When computing the realization shown in Figs. 1-5, we took iV mix = 1 and N v = 40. 
Because of the significant thermal broadening of the profile and the decrease in small-scale 
pulsation velocity amplitude with decreasing scale, the difference between the realizations 
for N v = 40 and N v ^> 40 is small (see the full paper for a discussion). The dependence 
on iV m i x is discussed in Section 7 below. 

We calculated the velocity field u(x, y, z) in a cube (2/) 3 with an edge 21. The function 
u(x, y, z) satisfies periodic boundary conditions at the boundaries of the large cube L 3 . In 
our calculations (Figs. 1-5), we took L = 1 and 21 = 0.4. We chose 21 < L to reduce the 
influence of the periodic boundary conditions. The cube (2/) 3 was covered by a 3D mesh 
of 160 3 computational points. The separations between the mesh points are 51 = 2//160. 
Each column in Fig. |21 has a length 21 along the line of sight and a cross section in the 
plane of the sky in the form of a 51 x 51 square. 

1.2 Turbulence and the Kinematic SZ Effect 

The hydrodynamic velocity distribution of a cluster in the plane of the sky can also be 
analyzed by studying manifestations of the kinematic SZ effect (Sunyaev and Zeldovich 
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Figure 7: Distribution of the CMB intensity fluctuations (1.4) in the y, z plane of the 
sky that follows from the Kolmogorov model, (a) A map of the line-of sight momentum 
p(y, z); L = 1. The white and black colors correspond to the motion away from and toward 
the observer, respectively. The shades of gray give the local values of AT r (y, z)/T r . (b) 
The AT r (y,z)/T r = isoline that separates the black and white regions with opposite 
signs of AT r . (c) The p(y, z = 0) cut (solid curve) of the p(y, z) surface. Dashed line 
represents the same cut in the case of the function p abs given by formula (1.6). 

1980), i.e. by measuring the intensity fluctuations of the cosmic microwave background 
radiation (CMB) within the cluster 



Sunyaev et al. (2003) and Nagai et al. (2003) used this method of analysis in their 
numerical modeling. Formula (1.4) was first suggested by Sunyaev and Zeldovich (1970) 
to calculate the primordial Doppler CMB fluctuations and was used by Sunyaev (1977) 
in calculating the fluctuations due to secondary ionization. In formula (1.4), y and z 
are the coordinates in the plane of the sky; x is the coordinate along the line of sight; 
and v\\ is the component of the local hydrodynamic velocity vector along the line of 
sight (the x velocity component) below denoted by u, v\\ = v\\/u rms ] and T r is the CMB 
temperature. The AT r /T r fluctuations (1.4) differ only by a factor from the fluctuations 
in the dimensionless total momentum 



of the emitting matter along the x axis (the line-of-sight or column momentum), where 
h e = n e /(n e ) centT , (n e ) centr , is the central cluster density. Below, the density n e is assumed 
to be roughly uniform over the cluster core. 

Fig. [7|and|Hlshow maps, reliefs, and cuts of the intensity fluctuations (1.4). The degree 
of blackening in Fig. [7| (a) is proportional to the local AT r (y, z)/T r fluctuation amplitude 
and the local momentum p(y,z). Fig. [7| (b) shows regions in the y, z plane of the sky 
that move 9 toward (black) and away from us (white). We performed our calculations by 
using the Kolmogorov model (for details, see the sections below). For long- wavelength 

9 Here, we are dealing with the motion on average, because p = J udx/(2l) 




(1.4) 




(1.5) 
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turbulence with iV m i x ~ 1 (see Fig. the size of the region with the same sign of p (1.5) 
is on the order of the cluster size L. In the typical realization shown in Figs. 1-5 and 
7-10, N mix = 1 and N„ = 40. Fig. [7| (c) shows the cut of the p(y, z) surface (see Fig. |HJ) 
by the z = plane that passes through the cluster center. The solid and dashed curves 
in this figure indicate, respectively, the functions p(y, 0) (1.5) and 

Pabs{y,0) = J n e \u(x,y,0)\.dx/(2l) (1.6) 

In the large-scale case (iV mix ~ 1), the functions \p(y,z)\ (1.5) and Pahs(y, z) (1.6) are of 
the same order of magnitude, as illustrated by Fig. [7| (c). 

The sizes of the identically signed spots (Fig. [7| (b) for the function p(y, z) and the 
\p\/Pahs ratio are determined by the scale L mix = L/N mix . If there are many mixers 
(-^mix 3> 1), then the function p(y,z) changes sign during a displacement Al± ~ L mix 
in the plane of the sky transverse to the line of sight. In this case, \p\/pabs ~ 1/V^mix- 
In deriving these estimates, we assume that (1) the separation between mixers is on the 
order of their size, (2) the fluctuation amplitude of the velocity u Tms on each mixer is of 
the same order of magnitude, and (3) the velocity correlation decays on a scale on the 
order of L mix . In this case p a b s ~ 1, because (u) ~ 1, (2/) ~ 1 (we set h e = 1), and 
\p\ ~ iV" mix , because the total momentum p(y, z) (1.5) along the x line of sight is the sum 
of the momenta of iV m i x out-of-phase mixers. 

It is interesting to calculate the Kolmogorov scaling law for the CMB fluctuations 
AT r /T r oc p (1.4). This law determines the amplitude of the fractal "fringes" (small 
fluctuations or spots) 10 in the dependences AT r (y, z) and p(y, z) (see Figs. [3(a), (c), and 
IHJ). The velocity is known (see Sections 3 and 4) to be given by 

(\v(r 2 ) - v(n)\) = (Av) Al oc (A/) 1/3 , Al = |r 2 - n|. (1.7) 

Let us show that the following law holds for the kinematic SZ effect: 

(|AT r (r ±2 ) - AT r (r xl )\) = [A(AT r )] Al± oc (Ap) Al± oc (Al ± f 6 , (1.8) 

where A/j_ = \r±2 — r±x\, r± = (y, z) is a two-dimensional vector in the plane of sky. We 
restrict our analysis to the case iV m i x ~ 1 presented in Figs. 1-5 and 7-10. Below, we 
give a brief derivation. Let us consider two adjacent parallel x lines of sight separated by 
a distance Alj_. The value of (Ap) A i ± is equal to the difference between the momenta p 
(1.5) on these lines of sight. Let us calculate this difference. The contributions of large 
eddies with a size Z rot (Z ro t > Al±) to the difference between integrals (1.5) are small, 
because these contributions are almost equal on adjacent lines of sight. The contribution 
of small eddies (Z rot < Al±) increases with size / rot and reaches a maximum at (Z rot ~ A/j_). 
Consequently, eddies with sizes on the order of the separation between the lines of sight 
should be considered to estimate (1.8). 

These eddies are uncorrelated. Indeed, the correlation in the velocity field induced by 
them decays on a scale ~ Al±. Therefore, the corresponding contribution to the difference 

10 Small oscillations or fringes are formed by a random addition of many eddies with different scales; 
see the sections below and the full paper. 
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Figure 8: Relief of the kinematic SZ effect. The height p of the point (y, z, p) on the p(y, z) 
surface (the height above the (y, z) plane) is proportional to the fluctuation amplitude 
AT r (y,z)/T r (lA). 



A(AT r ) oc Ap is accumulated through a random addition. According to formulas (1.5) 
and (1.7), the momentum from a single eddy is 

~ u Tms (Al ± /L mix )' 1/3 {AljL)/u rms . 

There are ~ L/Al± such eddies on a line of sight which has a length ~ L ~ ^mix- 
Expression (1.8) follows from this and from the fact that the effect is proportional to the 
square root of number of small-scale eddies under consideration, because 

(Ap) Al± ~ (AZ ± /£mi*)- 1/3 (AZ ± /L) ^L/AT L . 

Scaling (1.8) is valid in the Kolmogorov range of scales L v < Al± < L mix . In observations 
of AT r /T r oc p or in direct numerical simulations, the small-scale limit is determined by 
the aperture resolution or the spatial mesh step, respectively. 

1.3 Comparison of the Two Methods for Studying Turbulence 

Let us compare the two effects under discussion, the Doppler broadening of X-ray lines and 
the kinematic SZ effect. The map of CMB intensity (1.4) carries information about the 
momentum of a unit column along the line of sight. This gives the velocity relative to the 
frame of reference associated with the CMB near the cluster 11 , i.e. its peculiar velocity. 

11 The velocity of the cluster as a whole can be determined by integrating the line profile over the plane 
of the sky within the cluster and finding its centroid. This operation will allow us to measure the sum of 
the cluster recession velocity due to the expansion of the Universe (redshift) and the line-of-sight peculiar 
cluster velocity in the frame of reference associated with the CMB. The kinematic SZ effect allows this 
peculiar velocity component to be determined. 
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The X-ray line profiles are more informative. In the simplest case of an isothermal plasma 
with a constant iron abundance, they provide information about the velocity distribution 
of matter relative to the velocity of the cluster as a whole (in particular, the velocity 
dispersion). Below, we discuss in detail several characteristic features of the spectrum 
and their relationship to the turbulent velocity field. 

To compare the two effects, Fig. El shows a map of the velocity distributions fn(AE) 
over the (y, z) cluster plane. The distributions /g(A£?) were computed at the points 
located at the center, vertices, and middles of the sides of the 21 x 21 square shown in Figs. 
E|and|Hl We took a frame of reference in which the total momentum J J p(y, z)dydz of the 
emitting and scattering matter is equal to zero. We are looking for the distributions 
in this frame. The scale of dimensionless velocities u was recalculated to energy shifts AE 
by using formula (1.3) at u rms = 270 km s -1 . The velocity dispersion changes from point 
to point in the plane of the sky, u rms = u Tms (y, z). For definiteness, in the normalization 
M rms = 270 km s _1 , we chose the dispersion at the center y — 0, z — 0. Note that the 
dispersion over the cluster as a whole appreciably 12 exceeds the local dispersion. The 
relief of the kinematic SZ effect is determined by the shift velocity distribution (average 
velocity of the column as a whole). We see from Figs. El and El that these shift velocities 
give a large contribution to the velocity dispersion in the cluster as a whole. 

The /ht(A.E') profiles of Fe XXV (Eq = 6.7keV) lines with turbulent and thermal 
broadenings are drawn in Fig. They were obtained from the distributions fn{AE) in 
Fig. El by folding with the Maxwellian distribution at kT = 3 keV. Thermal smoothing 
"blurs" some of the features in the hydrodynamic distribution fn 13 , although the spectra 
remain fairly complex and varied. The centers of the fn{AE) and fm(AE) profiles 
coincide. Of course, they coincide with the profile center determined from the line-of- 
sight momentum p(y, z). Therefore, the map of the kinematic SZ effect (Figs. |7|and|Sl) is 
simultaneously the map of the shifts in the profile centers in Figs. El an dEIl Let us show 
this by using an example. The black color in the left panel of Fig. [7| (a) (the lower left 
corner of the square) corresponds to high velocities directed toward the observer. These 
velocities clearly show up in the rightward shift of the profile along the energy axis from 
AE = in the spectrum in the lower left corner of Figs. El and EH because motion toward 
the observer causes the photon energy to change. The shifts of the profile center as the 
position varies in the plane of the sky are large (3-6 local velocity dispersions). Thus, it 
is possible to independently and experimentally determine the shifts of the profile center 
by two different methods (the kinematic SZ effect and analysis of the line profile). 

The line shape completely changes during the shift to a large distance Al± ~ L in a 
direction transverse to the x line of sight (N mix ~ 1). The shape changes continuously: the 
change in shape is small for small shifts. For example, the variation of such a parameter 
as the shift of the profile center during transverse shifts is shown in Fig. [7| (c). The 
profiles significantly deviate from the Gaussian profiles (a large shift of the spectrum 
center, splitting of the maximum, a sharp break in the exponential tails). A comparison 
with the Gaussian distribution for the middle spectrum from the upper triplet of spectra 
(Figs. l2land irU|) is shown in Fig. [TJ 

The presence of turbulence in the intergalactic cluster gas and its amplitude can also 
be judged from differences in the intensity distributions of permitted from one side and 

12 by a factor of 1.5 to 2.5; cf. parabolas 2 and 4 in Fig. 
13 At u rms = 270 km s" 1 and kT = 3 keV, A^dh/A^dt = 2.7. 
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Figure 9: A map of the fn{AE, y, z) profiles for iron ions in the (y, z) plane at the points 
that coincide with the center, vertices, and middles of the sides of the square shown in 
Fig. 0(a). 
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Figure 10: The /ht(A£?) profiles of Fe XXV lines at kT = 3 keV (the same as Fig. El but 
in linear coordinates). The spectrum shown in the middle panel of the upper row from 
the three panels was used to construct Figs. [T] and |U (a). The central and bottom left 
panels were used in Figs. IU(b) and (c). 
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forbidden, intercombination lines and their satellites from another side over the cluster 
(Gilfanov et al. 1987). The optical depth in permitted lines can appreciably exceed unity 
in the absence of turbulence. Turbulent broadening reduces the optical depth in lines 
and this effect. This is a third independent method for studying turbulence in clusters of 
galaxies (Churazov et al., 2003). 

Note that Fabian et al. (2003) pointed to the existence of long cold if a -emitting 
filaments in the Perseus cluster of galaxies as an argument against the well-developed 
turbulence in this cluster. 

The full paper is arranged as follows. Section 2 gives a geometrical formulation of the 
problem with projection onto the observer's direction (line of sight). Section 3 presents 
three main components of the Kolmogorov model: two boundaries with the sizes (I/ m ix, 
L u ) of the inertial hierarchy of eddies and fluctuation amplitude scaling. Section 4 is 
devoted to the Fourier expansion with random phases of the turbulent velocity field. 
The procedure for calculating the velocity distribution fn{ u ) from the specific realization 
u(x) of a random velocity field is described in Section 5. The distribution fuiu) can be 
approximately calculated from the first moments of the turbulent velocity field (Section 
6). When there are many mixers, iV" mix ^> 1, the turbulent distribution tends to a 

Gaussian distribution f G (u) (Section 7). The superposition of thermal and hydrodynamic 
broadenings is analyzed in Section 8. A theory of the turbulent profile shape is constructed 
in Sections 9 and 10. We elucidate the question of typical spikes in the profile and its 
wings and the variety of possible shapes. The transverse correlations (in the plane of 
the sky) are studied in Section 11. The Lorentz line broadening is taken into account in 
Sections 12 and 13. 

2 Turbulence and Doppler Shift 

We are interested in the emission from the hot intergalactic plasma that fills cluster of 
galaxies. The plasma is in hydrostatic equilibrium in the gravitational well produced by 
collision-less dark matter. The density decreases toward the periphery of the well on scales 
of the order of the cluster scale L. The kinetic energy of the turbulent motion caused by 
mergers reaches ~ 15% of the thermal energy. By the meaning of the formulated problem 
(Doppler shift), the hydrodynamic density and plasma temperature fluctuations, which 
affect the local intensity but not the frequency, are of little importance. The velocity 
fluctuations are dominant. 

There is a three-dimensional velocity field v(f,t). It is necessary to determine the 
Doppler shift in a frame of reference in which the cluster as a whole is at rest. We will 
consider the line profile when the X-ray telescope is pointed at a point in the cluster plane 
of the sky. The profile is defined by the function u(x, t) that specifies the velocity u along 
the x line of sight passing through the observer and the point at which the telescope is 
pointed. The small eddies that are localized at the points separated along the x axis 
are statistically equivalent. The turnover time of large-scale eddies is much longer than 
the time it takes for light to pass through the cluster. In addition, the evolution time of 
even the smallest observed velocity fluctuations is much longer than the exposure time 
during which an object is observed from a satellite. Consequently, we may omit the time 
dependence of the function u(x,t) and deal with the instantaneous velocity field u(x). 
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The Doppler change in photon energy AE = E — E is 

AE/E = —u/c, AE [eV] = 6.7 (w/300 [kms" 1 ] ) (_E /6700 [eV] ) (2.1) 

where E is the energy of the photon emitted from a region that is at rest relative to the 
cluster as a whole, and E is the energy of the photon emerged from a region that moves 
with velocity u relative to the adopted frame of reference. Let the cluster be to the right 
of the observer on the x axis. Accordingly, the velocity u is positive if the emitting region 
recedes from the observer. 

3 The Kolmogorov Model 

An individual line profile corresponds to each velocity field u(x). We are interested in 
statistically representative realizations of u(x). The Kolmogorov views of turbulence are 
of great importance in constructing a typical profile. The point is that they solve the 
difficult question regarding the contribution of small-scale fluctuations. 

In a Kolmogorov cascade, energy is transferred from large to small scales. The largest 
eddies have scales on the order of the mixer size L m[x . The size L v of the smallest eddies 
is determined by viscous dissipation. A cluster of volume L 3 is assumed to be covered by 
a three-dimensional mesh of A^ ix statistically equivalent mixers, 

N mix = L/L mix . (3.1) 

The v(r) fluctuations are isotropic in velocity vector orientation. The x axis is not high- 
lighted in any way. Therefore, the Kolmogorov scaling (1.7) can be written for the x 
velocity component (the projection of three-dimensional fluctuations onto the line of 
sight). In the hierarchy of scales, the mean energy losses E within the inertial range 
are scale-independent (Kolmogorov 1941), 

E ~ [ ( Au) At f/At ~ [ (Au) Ax ] 3 /Ax, (3.2) 

where E are the hydrodynamic energy losses per gram of matter, and Au is the velocity 
variation due to the shift in time At or along the line of sight Ax. The Kolmogorov 
scaling follows from (3.2): 

(A^^-E^iAx) 1 / 3 . (3.3) 
The locally viscous dissipation is given by the formula 

E = (v/2)J2J2( dvi / dx i + 9vj/9xi) 2 ~ v [ {Au) v f/[ (Ax) u ] 2 , (3.4) 

where (Au) u and (Ai)„ are, respectively, the velocity and diameter scales of eddies on 
which viscous dissipation takes place. To estimate these scales, law (3.3) is extended to 
the range of dissipation through viscosity (Kolmogorov 1941; Monin and Yanglom 1965; 
Landau and Lifshitz 1986). We then find from formula (3.4) that 

(Ax), ~ u 3 ^/E^ ~ (c s /u rms f* Zf (3.5) 
N v = L/(Ax) u , (At) u ~ yfj/i, (3.6) 
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where c s is the speed of sound, u rms is the characteristic hydrodynamic velocity (rms 
deviation or dispersion) defined by the second moment 



(«rms) 2 = J U 2 dx/L 



of the function u(x) (the integral is taken over a range on the order of L), and l p is the 
proton mean free path. Large-scale eddies mainly contribute to the dispersion u rms and 
the kinetic energy of the turbulence. N u (3.5) gives the number of the smallest eddies on 
the characteristic cluster scale. Formulas (3.5) and (3.6) define the smallest scales. 

The size of small eddies (3.5) is mainly determined by the proton mean free path 
l p . The latter varies over a wide range (many orders of magnitude) with amplitude of 
the random magnetic field in the intergalactic plasma of the cluster. Below, we give the 
corresponding estimates. In the absence of a magnetic field, the plasma viscosity 

v = 0.48 v Tp l P = 1.5 • 10 27 (kT) 5/2 /n_ 2 , 



l P = v Tp r p = 7.2 • 10 19 (kT) 2 /n. 2 , t>Tp = ^2kT/m p 

is determined by Coulomb collisions between ions (see, e.g., Rosenbluth and Sagdeev 
1983); in what follows, the kinematic viscosity v and the mean free path l p are in CGS 
units, the temperature kT is in keV, n = 10 _2 n_ 2 cm -3 , m p is the proton mass. The 
gyroviscosity of plasma with a magnetic field is given by the expression (Rosenbluth and 
Sagdeev 1983) 

v = 0.5 kT/m p u cp = 0.25fT P ^L, tl = m p VT P c/eH, 

where u cp is the cyclotron frequency of protons, and is their Larmor radius. The 
viscous scale ratio is 

|g|g j= jjj = (1-92 r p co cp f* = 2.3 • 10 6 ( H_ 8 ) 3 / 4 ( kT ) 9 / 8 ( n_ 2 )" 3 / 4 , 

where H = 10~ 8 H_ 8 G. Accordingly, the dimensionless viscous wave numbers N v (3.5) 
are 

Nu {H = 0) = 900 ( AUx ) 1/4 Ma 3 / 4 ( L 200 ) 3/4 ( n„ 2 ) 3 / 4 ( kT)~ 3 ^ 2 , 

N V {H + 0) = 3.5 • 10 9 ( N mix ) x / 4 Ma 3 / 4 ( L 200 ) 3/4 ( H. % ) 3/4 ( kT )" 3 / 8 , 

where Ma = u rms /c s ; u Tias and c s are the characteristic hydrodynamic velocity and the 
speed of sound, respectively; and L = 200L 200 kpc. 

In a hot rarefied plasma without a magnetic field, the mean free path l p is very large. 
Thus, a situation where N v is very small (iVj, ~ 1) is hypothetically possible. As we 
show in Sections 9 and 10 of the full paper, the line profile in this situation exhibits 
features that allow it to be distinguished from the case with N u ^> 1. Currently available 
estimates based on Faraday rotation measurements give magnetic field strengths in the 
range 1-10 /xG (see, e.g., Ge and Owen 1993), suggesting that the Kolmogorov range in 
the intergalactic turbulence spectrum is very wide. 
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4 Fourier Expansion 



For the line profile to be constructed, we must know the velocity distribution fn{u). The 
distribution fuiu) can be calculated from the velocity field u(x). Let us write out an 
analytical formula for the function u(x) that satisfies law (3.3). It would be natural to 
represent the stochastic turbulent motions in the Kolmogorov range of scales as a Fourier 
series 

N v 

u'(x; f) = ( N mix T Re n~ a ~ 1/2 exp(iip n ) exp (i2irnx/L) . (4.1) 

In formula (4.1), the index a specifies the velocity fluctuation scaling (Au)ax oc (Ax) a . 

— * 

In case (3.3), a = 1/3. The phases ip = {ip n } in the harmonic expansion (4.1) are 
independent random variables uniformly distributed in the segment [0,27r]. The velocity 
is given by the real part of expansion (4.1) (the symbol Re). The summation limits (4.1) 
coincide with boundaries (3.1) and (3.5) of the inertial range. The function of x (4.1) is 
considered in the segment [I1J2] with a length on the order of L. Series (4.1) describes 
the fractal curve u(x) whose fluctuations satisfy the Kolmogorov scaling law (3.3). This 
series makes it possible to carry out specific calculations. 

Let us "displace" the auxiliary function u' (4.1) so that the velocity field is centered 
(/ Mo d x = 0). We change u' — > uq, where 

rfo 

u (x) —u'(x) —ui, ui — I u'(x) dx/ki, ki — k — h, (h,k) ~ L, (4.2) 

Jh 

the function u' is calculated from formula (4.1). We normalize the function uq (4.2) to the 
rms velocity u Tias . By definition, u rms gives the dispersion a of the velocity distribution. 
The velocity dispersion, along with the kinetic energy of the turbulence, is determined by 
large eddies. We write 

u{x) = u (x)/u ims , M r 2 ms = / [u (x)] 2 dx/l 2 i. (4.3) 

Jh 

Below, we deal with the centered (4.2) and normalized (4.3) function u(x) defined by 
formulas (4.1)-(4.3). The model of a homogeneous cloud with size L is used as the first 
approximation to the density distribution. Clearly, this picture can be easily generalized 
to the real density distribution in clusters of galaxies. 



5 The Velocity Distribution 

Let us explain the meaning of the distribution /n(tt). There is a function u(x). Let us 
divide the u axis into segments 5u that are small compared to the viscous fluctuation 
amplitude (Au) u . Consider an arbitrary point u± on the u axis. By definition, 

/#(«*) = dp/du = ^^5x/l 2 iSu, (5.1) 

where the differential dp characterizes the "weight" of the set of subregions that move 
with velocity u*. We assume that the argument of the function u(x) belongs to the region 
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h < x < l 2 inside the cluster. The condition u± — 5u/2 < u(x+) < + 5u/2 is satisfied 
in each of the small segments 5x in sum (5.1). In this condition, the points x* belong 
to any of the small segments dx in sum (5.1). Because of the statistical pulsations in 
velocity u(x), sum (5.1) can contain many terms. Their number is approximately equal 
to the numbers of intersections of the function u = u(x) with the straight line u = u±. 
It is easy to see that distribution (5.1) is normalized, fn du = 1. The corresponding 
illustration is presented in Fig. [TT] 

Let us first consider an example with a monotonic function u(x) and then return to 
the non-monotonicity and multivaluedness of the function x(u) that is the inverse of u(x) 
and to the turning points. In the monotonic case, 



l/[du(x)/dx] = dx{u)/du = hifiiiu); 



7 du(x) 

hi— — 



dx 



i.e., the velocity gradient determines fn- The equations u = u(x) and fn 

parametrically specify the distribution /#(«)• Eliminating the parameter x from this pair 
of equations yields an explicit dependence fuiu). 

Let u{x) = U sin(27rx/L). Eliminating the parameter x, we then obtain 



L 1 



f H (u) = (hidu/dx) =——- — t-- = — = . (5.2) 

2nl 21 U cos(2nx/L) 2nl 2l ^/U 2 - u 2 

Expression (5.2) specifies the dependence in explicit form. Of course, an explicit 

expression can be derived only for very simple functions u(x). 

The example with an oscillating function u(x) clearly illustrates an interesting singu- 
larity of the distribution fu in the cold 14 case. It becomes infinite 15 at the stop, turning, 
or cuspidal points of the inverse function x(u) that are the extrema of the function u(x); 
because of the root behavior, the distribution /# (w) is asymmetric about the singularity. 
The case with one harmonic (5.2) is an example that greatly differs from the Gaussian 
distribution. In this example, the dependence Jh{u) has a minimum at the center and 
a singular break in the distribution tail. Consequently, in this example, the regions in 
which there is motion relative to the observer are more representative than the regions 
at rest (that is why the distribution maximum is off the center). As we will see below, 
in the case of many harmonics, the turning points or extrema lead to a needle-shaped 
distribution fniu). 

Thus, it is clear that distribution (5.1) is, in a sense, obtained by "projecting" the 
inverse function x(u) onto the u axis, as shown in Fig. ^2 (the connection between Figs, 
im (a) and (b)). Fig. ^2 (a) shows a realization of u(x) (4.1)-(4.3) with random phases 
ipn (N <n< N u , N mix = 2,N U = 100) in the interval h < x < l 2 {h = 0.2L, l 2 = 0.7L). 
The viscous limit is fairly close: (Ax) u = l 2 \/N u ~ 10~ 2 L, (Au) u ~ 0.1w rms . The viscous 
values determine the heights and widths of the u(x) spikes (Fig. II 

In calculating the distribution /#(«), we substituted the velocity u(x) with the set 
of its values Ui at points X{ = iSx (discretization). The Sx step was chosen to be small 
compared to the viscous length (Ax)j, in order to resolve dissipating eddies. The distri- 
bution fn(u) was calculated numerically. It was replaced by a histogram (Fig. [TT] (b)) 

1 In Section 8 of the full paper, we show how thermal Doppler broadening "smears" the singularities 
caused by turning points. 
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Of course, this is an integrable singularity. 
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Figure 11: A scheme for recalculating the instantaneous velocity field u(x) (a) to the 
velocity distribution fn(u) (b). The recalculation is performed by "projecting" the multi- 
valued function x(u) (the inverse of the function u(x)) onto the u axis. The fuiu) profile 
is similar to the profile of a hydrodynamically (without thermal broadening) broadened 
line. 
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with a small velocity step Su <C (Aw)„ in order not to lose any fluctuation. Because of 
the discretization, the above singularities of the distribution are cut off at turning 

points. The height of the spikes in the distribution is on the order of a/ (Am) v /5m,; 

cf. (5.2). The value of /(w*) at the step 

u± — 5u/2 < u < u* + 8u/2 

of the histogram is proportional to the number of discretization points Ui at the limits 
± Su/2. After the normalization J fudu = 1, we obtain the distribution /# shown in 
Fig.EKb). 

It would be natural to analyze the deviations from the Gaussian distribution together 
with the calculation of moments, because the extent to which the distribution is non- 
Gaussian can be characterized by the deviations of moments from their Gaussian values. 



6 The Method of Moments 

Above, we described the procedure for an exact calculation of the distribution f H (the ex- 
ample in Fig. ITTj) . Concurrently, it is instructive to describe the method for its analytical 
calculation. We approximate the distribution fjj by a polynomial 

f(u) = a + diu + a 2 u 2 + ... = djV? . (6.1) 

3=0 

The calculation procedure is based on the fact that the moments are determined by two 
independent methods. Let us calculate the moments \i n directly from the velocity field 
u(x) and the distribution fuiu). By definition, 

pl% POO 

fi n = [u(x)] n dx/l 21 = / u n f{u)du. (6.2) 



The first three moments are known (do not depend on the phases ip n ): fiQ = 1 (the 
normalization condition f f H du = 1), /ii = (centering (4.2)), and // 2 = 1 (the choice of 
velocity unit (4.3)). 

We truncate series (6.1) at the approximation order n app . Expansion (6.1) contains 
n &x>p + 1 unknown coefficients a n . Let us set up a system of equations to determine them. 
The equations are fi™ = \d n , where the expressions for and \d n are calculated from u and 
fn, respectively (see (6.2)). Let us first calculate the missing moments /Xg,/^,, fi^-app 
from the velocity u(x). For a given realization, these are just some specific numbers. We 
then express integrals (6.2), 

ti= ru n f(u)du, (6.3) 



in terms of expansion (6.1). 

Moments (6.3) are linear forms in unknowns a n . The coefficients of the forms depend 
on the left, u\ < 0, and right, u r > 0, integration limits in (6.3). Clearly, the power-law 
approximation (6.1) does not describe the decaying tails of the distribution f H . Polyno- 
mial (6.1) intersects the u axis at points u\ and u r . Inside the segment ui < u < u r , the 
polynomial is positive. 
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Figure 12: Comparison of the distributions fuiu) : the thin fluctuating solid curve - 
exact calculation, the thick smooth solid curve - calculation by the methods of moments 
(a polynomial up to u b inclusive), the thick dot-dashed curve - a polynomial up to u 6 , the 
thin dashed curve - Gaussian distribution f G . The distribution f G is symmetric about 
the vertical line u — 0. 

Solving the linear system of equations \i u n = \J n for the unknowns a n , we express a n in 
terms of the zeros ui and u T of the polynomial fn in the form of a rational 16 function of ui 
and u T . Substituting the derived expressions for a n (ui, u r ) into (6.1) yields the distribution 
fH{ui,u r ;u). The zeros ui and u r are defined by a system of two equations fH(ui,u r ;ui) 
= and fH(ui,u r ;u r ) = 0. At n app = 2 (approximation (6.1) is a parabola), ui = —V5 
and u r = v5- In high orders, the equations rapidly become cumbersome. 

Examples of polynomials for n app = 5 and 6 are given in Fig. El They refer to 
the realization of u(x) shown in Fig. ^2 The polynomials satisfactorily fit the exact 
distribution fn- The sixth-order polynomial (heavy dash-dotted line) catches even the 
additional peak in the right wing. Of course, polynomials of a limited order smooth out 
the "needles" of the exact distribution. The Gaussian distribution is indicated in Fig. 
1721 by the dashed curve. We see the essentially non-Gaussian behavior, which manifest 
itself in the profile asymmetry relative to the dashed u = vertical line (Fig. IT2j) and in 
the break of the distribution tails and the appearance of additional peaks. In addition, 

16 The ratio of polynomials 
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the odd moments are nonzero, and the ratios of sequential even moments differ markedly 
from the Gaussian ratios. For the typical realization shown in Fig. ^2 A*3 — 0.51 and /X5 
= 2.97. The measure of smallness of the odd moments is determined by the normalization 
(j,2 — 1 (4.3). The increase in the absolute value of the moment with increasing moment 
number should also be taken into account. In the Gaussian case, the following relations 
hold: /if = 3 (/if) 2 , = 15 (/if) 3 . In our case, however, /X4/3 = 0.93 and fie/15 = 0.71. 

Let us consider the influence of viscosity. To assess the role of the width of the 
inertial range logio^-, we compared the velocities u(x) and the distributions fuiu) with 
N v ~ 10 2 and N v ~ 10 3 . As should be the case, the differences are attributable to the 
appearance of small-scale fluctuations as N v (3.5) greatly increases. The positions and 
shapes of the significant fluctuations with A u ~ 0.lM rms depend weakly on N u if N v > 100. 
The first moments (up to the sixth moment inclusive) change only slightly (by ~ 1%). 



7 Gaussian Asymptotics 

The main regulator on which the deviation of /# from the Gaussian distribution f G 
depends is the parameter N m[x that specifies the extent to which the scale of turbulence 
(3.1) is large. The viscous scale in this ratio is unimportant. Therefore, we set A^ = 00. 
The difference between the distributions and f G decreases with increasing N mix . At 
A^mix — 00, the profile is Gaussian. Let us first show this and then present the results 
that illustrate the pattern of convergence fn — » f G with increasing N m ^. 

Let us analyze the case with N m i x ^> 1. Above, the calculations were associated with 
the dependence of velocity u on coordinate x. In this case, the phases ipn (4.1) were fixed. 
This is an approach with averaging over the coordinate. Let us look at the question 
differently. Let x be fixed (x = xg x ). The quantity uix^iji) = u(ip) gives the velocity 
at point xrfix for a given set of phases ip (4.1). Let us study the statistics of u{ip) (an 
approach with phase averaging). 

Since the phases ip n are independent and random and the position x is fixed, we may 
omit the regular oscillating factors exp (i2iinx/L) in (4.1). Consequently, the velocity 
u(x &x ; V) is given by 

00 

« = (iV mix ) 1 /3 Un 5/6 , (7.1) 

™=A f mix 

where, in place of the phases ip n we introduced the equivalent independent random vari- 
ables £ n (£ n ) = 0, and (£ 2 ) = 1; () denotes averaging. We assume, for simplicity, that ip n 
takes on values of +1 or -1 with equal probabilities. The probability distribution function 
of the random variable u (7.1) should be determined. 
We see from definition (7.1) that the mean velocity (u) is equal to zero. The probability 
distribution /h(w) is symmetric (/# is an even function of u). All odd moments are zero 17 . 
It is easy to show that the second moment of (7.1) is 

00 

(n 2 ) = (iV mix ) 2 / 3 ^ 5/3 = (A r mix) 2/3 C(5/3,iV mix ), (7.2) 

n=AT mix 

17 In the specific realizations considered in Sections 4-6, the function is asymmetric. Here, we 

deal with the probability distribution after averaging over all possible phases (i.e. over all possible 
realizations). This distribution is symmetric. 
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where C(5/3, iV m i x ) is the generalized Riemann zeta function (Vinogradov 1979). In the 
limit iVmix > 1, (u 2 ) = 3/2, because ^n~ 5//3 « J dn/n 5 ^ 3 . 
Let us calculate the fourth moment. We have 

«V( N mix ) 4 / 3 = ( J2 Un 5/6 ) 4 = -2j2 ^/n w/s + 3 ( £ £> 5 / 3 ) 2 + A (7.3) 

In (7.3), the term A is the sum of the terms with an odd number of factors £ n . After 
averaging, the term A vanishes, and we obtain 

(« 4 )/( N mix ) 4 / 3 = -2 n ~ W/3 + 3 ( E ^' 5/3 ) 2 = -2 C(10/3, iV mix ) + 3 [ C(5/3, iV mix ) ] 2 , 

(7.4) 

because £ n = ±1 and £ 4 = 1. For A^" mix ^> 1, expression (7.4) tends to the limit 

<« 4 ) = 3 (3/2) 2 . (7.5) 

As we see, the ratio of the fourth moment (7.5) to the square of the second moment (7.2), 
(it 4 )/( (u 2 ) ) 2 , is equal to three. This should be the case for the Gaussian distribution. 
The aforesaid is also true for the random variables £ n of a more general form than £ n = 
±1. Indeed, for A^" mix ^> 1, sum (7.1) includes many approximately equal (in absolute 
value) independent random terms. In this case, in view of the central limit theorem, the 
probability distribution of the random variable u is Gaussian. At A^ mix ~ 1, this is not 
the case, because th first terms in sum (7.1) significantly differ due to the factor n -5 / 6 . 
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-0.16 
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It remains to study the convergence f H F with increasing iV m i x . The corresponding 
results are presented in Table 2. 

This table shows how the moments converge to their Gaussian values (the last row). 
The moments fi n are defined by formulas (6.2), R4 = /14/3/i 2 ,, R& = He/l-Sfi 3 ,- The first 
moment is equal to zero due to centering (4.2). Because of normalization (4.3), the second 
moment is always equal to unity. In the Gaussian limit, fi3 = ^5 = and R4 — R% — 1. 
The procedure for calculating the moments /i 3 — /i 6 was described in Section 5 and 6. The 
table list the mean values obtained by averaging over many realizations l;,u(x) (4.1). For 
the odd moments, we average the absolute value of the moment to establish the degree 
of deviation from zero. The realization-averaged values are given in rows 3-5 of Table 2. 
For comparison with a unit realization, the second row gives the values that refer to the 
example in Fig. ^2 see a ls° the end of Section 6. This row is highlighted by the symbol 
(*). The formula 

Rl(N mix ) - 1 = (n 4 )/3( (u 2 ) ) 2 - 1 = (-2/3) C(10/3, N mix )/[ C(5/3, N mix ) ] 2 , 
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used in column 5 of the table follows from expressions (7.2) and (7.4). It shows the 
deviation from the Gaussian distribution in the case of averaging over all realization £ n . 
The corresponding i? 4 calculated form the zeta function is denoted by the superscript (. 

Let us discuss the data in the Table 2. We see that the moments approach their 
Gaussian values as iV mix increases (refinement of the leading or dominant turbulence 
scale (3.1)). Therefor, the distribution f H also tends to f G . As should be the case, the 
convergence in lower moments (^3,^4) is faster than the convergence in higher moments 
(fj,5,fj,e). The deviations in higher moments are larger, and they decrease with increasing 
iV m i x more slowly. The main conclusion is that for large-scale turbulence with A^ m ; x = 2—10 
{L m i x ~ L, (3.1)), there are significant deviations from the Gaussian distribution. 
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